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Abstract: We discuss the basic difficulties in understanding the origin of the highest energy particles in 
the Universe - the ultrahigh energy cosmic rays (UHECR). It is difficult to imagine the sources they are 
accelerated in. Because of the strong attenuation of UHECR on their propagation from the sources to us 
these sources should be at cosmologically short distance from us but are currently not identified. We also 
give information of the most recent experimental results including the ones reported at this conference 
and compare them to models of the UHECR origin. 



Introduction 

More than forty years ago, in 1963, John Lins- 
ley [1] published an article about the detection of a 
cosmic ray of energy 10^" eV. The article did not 
go unnoticed, neither it provoked many comments. 
The few physicists that were interested in high en- 
ergy cosmic rays were at that time convinced that 
the cosmic ray energy spectrum will continue for- 
, ever. The fact that cosmic rays may have energies 
' exceeding 10^ GeV (10^'^ eV) was established in 
I the late thirties by Pierre Auger and his collabora- 
' tors. Showers of higher and higher energy were de- 
, tected in the mean time - seeing a 10^° eV shower 
' seem to be only a matter of time and exposure. Al- 
' ready in the fifties there was a discussion about the 
origin of such ultra high energy cosmic rays and 
I Cocconi [2] reached the conclusion that they must 
be of extragalactic origin since the galactic mag- 
netic fields are not strong enough to contain such 
particles. 

How exclusive this event is became obvious three 
years later, after the discovery of the microwave 
background radiation (MBR). Almost simultane- 
ously Greisen in US [3] and Zatsepin&Kuzmin [4] 
in the USSR published papers discussing the prop- 
agation of ultra high energy particles in extragalac- 
tic space. They calculated the energy loss dis- 
tance of nucleons interacting in the microwave 
background and reached the conclusion that it is 



shorter than the distances between powerful galax- 
ies. The cosmic ray spectrum should thus have 
an end around energy of 5x10^^ eV. This effect 
is now known as the GZK cutoff. 

The experimental statistics of such events grew 
with the years, although not very fast. The flux of 
UHECR of energy above 10^° eV was estimated 
to be 0.5 to 1 event per square kilometer per cen- 
tury per steradian. Even big detectors of area tens 
of km^ would only detect few events for ten years 
of work. The topic became one of common inter- 
est during the last decade of the last century when 
ideas appeared for construction of detectors with 
effective areas in thousands of km'^ [5, 6]. Such 
detectors would detect tens of events per year and 
finally solve all mysteries surrounding UHECR. 
The Auger observatory (3,000 km^) is now almost 
fully deployed and the Telescope Array (TA, 1,000 
km^) is being constructed in Utah, U.S.A. The ex- 
pectations for the flux of UHECR are now smaller. 
The High Resolution Fly's Eye (HiRes) and Auger 
have shown that the rate of events above 10^" eV 
is about 10 times smaller than previously thought. 

Cosmic rays are defined as charged nuclei that 
originate outside the solar system. They come on 
a featureless, power law like, F{E) = K x E'", 
spectrum that extends beyond 10^^ GeV per par- 
ticle. There are only two distinct features in the 
whole spectrum. At energy above 10^ GeV the 
power law index a steepens from 2.7 to about 3.1. 
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This is called the knee of the cosmic ray spectrum. 
At energy above 3x10^ GeV the spectrum flattens 
again at the ankle. 

The common wisdom is that cosmic rays below the 
knee are accelerated at galactic supernovae rem- 
nants. Cosmic rays above the knee are also thought 
to be of galactic origin, although there is no clue of 
their acceleration sites. Cosmic rays above the an- 
kle are thought to be extragalactic. 

More recently, with the improved accuracy and ex- 
posure of the modern experiments, several theoret- 
ical models that fit the measured ultra high energy 
cosmic rays (UHECR) spectrum have been devel- 
oped. We will discuss the experimental data and 
compare them to some of these models. 

Air Shower Detection Methods. 

Cosmic rays of energy above 10^^ eV are detected 

by the showers they generate in the atmosphere. 
The atmosphere contains more than ten interac- 
tion lengths even in vertical direction and is much 
deeper for particles that enter it under higher zenith 
angles. It is thus a deep calorimeter in which the 
showers develop, reach their maximum, and then 
start being absorbed. There are generally two types 
of air shower detectors: air shower arrays and opti- 
cal detectors. Air shower arrays consist of numer- 
ous particle detectors that cover large area. The 
shower triggers the array by coincidental hits in 
many detectors. The most numerous particles in 
an air shower are electrons, positrons and photons. 
The shower also contains muons, that are about 
10% of all shower particles, and hadrons. 

The direction of the primary particle can be recon- 
structed quite well from the timing of the differ- 
ent hits, but the shower energy requires extensive 
Monte Carlo work with hadronic interaction mod- 
els that are extended orders of magnitude above 
the accelerator energy range. The main composi- 
tion sensitive variable is the ratio of the number of 
muons in the shower to the number of electrons 
Ne, or the ratio of muon to electron densities at a 
certain distance from the shower axis. The type of 
the primary particles can only be studied in statis- 
tically big samples because of the fluctuations in 
the individual shower development. Even then it is 
strongly affected by the differences in the hadronic 
interaction models. 



The optical method uses the fact that part of the 
particle ionization loss is in the form of visible 
light. AH charged particles emit in air Cherenkov 
light in a narrow cone around their direction. In 
addition to that charged particles excite Nitro- 
gen atoms in the atmosphere, which emit fluo- 
rescence light. The output is not large, about 4 
photons per meter, but the number of shower par- 
ticles in UHECR showers is very large, and the 
shower can be seen from distances exceeding 30 
km. The fluorescence detection is very suitable 
for UHECR showers because the light is emitted 
isotropically and can be detected independently of 
the shower direction. Since optical detectors fol- 
low the shower track, the direction of the primary 
cosmic ray is also relatively easy. The energy 
of the primary particles is deduced from the total 
number of particles in the shower development or 
from the number of particles at shower maximum. 
The rough number is that every particle at maxi- 
mum carries about 1 .5 GeV of primary energy. The 
mass of the primary cosmic ray nucleus is studied 
by the depth of shower maximum X^axt which is 
proportional to the logarithm of the primary energy 
per nucleon. 

The Highest Energy Cosmic Ray Event 
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Figure 1: The shower profile of the highest energy 
cosmic ray shower detected by the Fly's Eye. 

The highest energy cosmic ray particle was de- 
tected by the Fly's Eye experiment [7]. We will 
briefly describe this event to give the reader an 
idea about these giant air showers. The energy of 
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this shower is estimated to be 3 x 10 eV. This is 
an enormous macroscopic energy. 3x10^" eV is 
equivalent to 4.8 x 10® erg, 7.2 x lO'^'' Hz or the en- 
ergy of 290 km/h tennis ball in the favorite units 
of Alan Watson. Fly's Eye was the first air fluo- 
rescence experiment, located in the state of Utah, 
U.S.A. Fig. 1 shows the shower profile of this event 
as measured by the Fly's Eye. Note that the max- 
imum of this shower contains more than 2x 10^^ 
electrons and positrons. Both the integral of this 
shower profile and the number of particles at max- 
imum give about the same energy. 

The errors of the estimates come from the errors of 
the individual data points, but mostly from the un- 
certainty in the distance between the detector and 
the shower axis. The minimum energy of about 
2x 10^° eV was calculated in the assumption that 
the shower axis was much closer to the detector 
than the data analysis derived. 

ORIGIN OF UHECR 

The first problem with the ultra high energy cos- 
mic rays is that it is very difficult to imagine what 
their origin is. We have a standard theory for the 
acceleration of cosmic rays of energy below the 
knee of the cosmic ray spectrum at galactic super- 
nova remnants. This suggestion was first made by 
Ginzburg&Syrovatskii in 1960's on the basis of en- 
ergetics. The estimate was that a small fraction (5- 
10%) of the kinetic energy of galactic supernova 
remnants is sufficient to maintain the energy car- 
ried by the galactic cosmic rays. The acceleration 
process was assumed to be stochastic, Fermi type, 
acceleration that was later replaced with the more 
efficient acceleration at astrophysical shocks. 

This statement stands, but it is not applicable to 
all cosmic rays. Much more exact recent estimates 
and calculations show that the maximum energy 
achievable in acceleration on supernova remnant 
shocks is not higher than 10^ GeV. This excludes 
not only UHECR, but also the higher energy galac- 
tic cosmic rays, that require supernova remnants 
in special environments [8]. There are now some 
very interesting ideas about shock magnetic fields 
amplification by cosmic rays that lead to higher ac- 
celeration energy and flatter energy spectrum. See 



the talk of Pasquale Blasi for a discussion of new 
acceleration models. 
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Figure 2: Requirements for acceleration of charged 
nuclei at astrophysical objects as spelled out by 
A.M. Hiflas (see text). 

The reader should note that currently the acceler- 
ation of charged nuclei at supernova remnants is 
mostly a theoretical prediction. Supernova rem- 
nants have higher matter density than interstellar 
space and one expects that the accelerated nuclei 
would interact with the matter and generate high 
energy 7-rays. 

Although Cherenkov gamma ray telescopes have 
observed supernova remnants with TeV 7-ray 
emission, there is no proof that TeV and higher en- 
ergy 7-rays are generated in hadronic interactions. 
On the other hand, multi-wavelength observations 
show the existence of very strong shocks in super- 
nova remnants with TeV 7-ray emission. 

We should then turn to extragalactic objects for ac- 
celeration to energies exceeding 10^° eV. The scale 
for such acceleration was set up by Hillas [9] from 
basic dimensional arguments. The first require- 
ment for acceleration of charged nuclei in any type 
of object is that the magnetic field of the object 
contains the accelerated nucleus within the object 
itself. One can thus calculate a maximum theo- 
retical acceleration energy, that does not include 
and efficiency factor, as Emax < jeZBR , where 
7 is the Lorentz factor of the shock matter, Z is 
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the charge of the nucleus, B is the magnetic field 
value, and R is the linear dimension of the object. 

Figure 2, which is a redrawn version of the original 
figure of Hillas, shows what are the requirements 
for acceleration to more than lO'^^ eV. The lower 
edge of the shaded area shows the minimum mag- 
netic field value for acceleration of iron nuclei as a 
function of the dimension of the astrophysical ob- 
ject. The upper edge is for acceleration of protons. 

There are very few objects that can, even before an 
account for efficiency, reach that energy: highly 
magnetized neutron stars, active galactic nuclei, 
lobes of giant radiogalaxies, and possibly Gpc size 
shocks from structure formation. Other potential 
acceleration sites, gamma ray bursts, are not in- 
cluded in the figure because of the time depen- 
dence of magnetic field and dimension. 

Possible Astrophysical Sources of UHECR 

In this subsection we give a brief description of 
some of the models for UHECR acceleration at 
specific astrophysical objects. For a more complete 
discussion one should consult a review paper on 
the astrophysical origin of UHECR [10], that con- 
tains an exhaustive list of references to particular 
models. 

• Pulsars: Young magnetized neutron stars 
with surface magnetic fields of 10^^ Gauss 
can accelerate charged iron nuclei up to en- 
ergies of 10^" eV [11]. The acceleration pro- 
cess is magnetohydrodynamic, rather than 
stochastic as it is at astrophysical shocks. 
The acceleration spectrum is very flat pro- 
portional to \/E. It is possible that a large 
fraction of the observed UHECR are acceler- 
ated in our own Galaxy. There are also mod- 
els for UHECR acceleration at magnetars, 
neutron stars with surface magnetic fields up 
to 10^^ Gauss. 

• Active Galactic Nuclei: As acceleration site 
of UHECR jets [12] of AGN have the advan- 
tage that acceleration on the jet frame could 
have maximum energy smaller than these of 
the observed UHECR by l/F, the Lorentz 
factor of the jet. The main problem with 
such models is most probably the adiabatic 



deceleration of the particles when the jet ve- 
locity starts slowing down. 

• Gamma Ray Bursts: GRBs are obviously 
the most energetic processes that we know 
about. The jet Lorentz factors needed to 
model the GRB emission are of order 100 
to 1000. These models became popular with 
the realization that the arrival directions of 
the two most energetic cosmic rays coin- 
cide with the error circles of two power- 
ful GRB. Different theories put the acceler- 
ation site at the inner [13] or the outer [14] 
GRB shock. To explain the observed UHE- 
CRs with GRBs one needs fairly high cur- 
rent GRB activity, while most of the GRB 
with determined redshifts are at Z > 1 . 

• Giant Radio Galaxies: One of the first con- 
crete model for UHECR acceleration is that 
of Rachen&Biermann, that dealt with accel- 
eration at FR II galaxies [15]. Cosmic rays 
are accelerated at the 'red spots', the termi- 
nation shocks of the jets that extend at more 
than 100 Kpc. The magnetic fields inside the 
red spots seem to be sufficient for accelera- 
tion up to lO^*' eV, and the fact that these 
shocks are already inside the extragalactic 
space and there will be no adiabatic decel- 
eration. Possible cosmologically nearby ob- 
jects include Cen A (distance of 5 Mpc) and 
M87 in the Virgo cluster (distance of 18 
Mpc). 

• Quiet Black holes: These are very mas- 
sive quiet black holes, remnants of quasars, 
as acceleration sites [16]. Such remnants 
could be located as close as 50 Mpc from our 
Galaxy. These objects are not active at radio 
frequencies, but, if massive enough, could 
do the job. Acceleration to lO^'' requires a 
mass of 10^ Mq. 

• CoUiding Galaxies: These systems are at- 
tractive with the numerous shocks and mag- 
netic fields of order 20 /uG that have been ob- 
served in them [17]. The sizes of the collid- 
ing galaxies are very different and with the 
observed high fields may exceed the gyrora- 
dius of the accelerated cosmic ray. 
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• Clusters of Galaxies: Magnetic fields of or- 
der several /zG have been observed at length- 
scales of 500 Kpc. Acceleration to almost 
10^° eV is possible, but most of the lower 
energy cosmic rays will be contained in the 
cluster forever and only the highest energy 
particles will be able to escape [18]. 

• Gpc scale shocks from structure forma- 
tion: A combination of Gpc scales with 1 
nG magnetic field satisfies the Hillas cri- 
terion, however the acceleration at such 
shocks could be much too slow, and subject 
to large energy loss. 

Top-Down Scenarios 

Since it became obvious that the astrophysical 
acceleration up to 10-^° eV and beyond is very 
difficult and unlikely, a large number of particle 
physics scenarios were discussed as explanations 
of the origin of UHECR. To distinguish them from 
the acceleration (bottom-up) processes they were 
called top-down. The basic idea is that very mas- 
sive (GUT scale) X particles decay and the result- 
ing fragmentation process downgrades the energy 
to generate the observed UHECR. Since the ob- 
served cosmic rays have energies orders of mag- 
nitude lower than the X particle mass, there are 
no problems with achieving the necessary energy 
scale. The energy content of UHECR is not very 
high, and the X particles do not have to be a large 
fraction of the dark matter There is a large number 
of topological defect models which are extensively 
reviewed in Ref. [19]. 

There are two distinct branches of such theo- 
ries. One of them involves the emission of X 
particles by topological defects. This type of 
models follows the early work of C.T. Hill [20] 
who described the emission from annihilat- 
ing monopole/antimonopole pair, which forms a 
metastable monopolonium. The emission of mas- 
sive X particles is possible by superconducting 
cosmic string loops as well as from cusp evapora- 
tion in normal cosmic strings and from intersect- 
ing cosmic strings. The X particles then decay 
in quarks and leptons. The quarks hadronize in 
baryons and mesons, that decay themselves along 
their decay chains. The end result is a number of 
nucleons, and much greater (about a factor of 30 in 



different hadronization models) and approximately 
equal number of 7-rays and neutrinos. 

A monopole is about 40 heavier than a X particle, 
so every monolonium can emit 80 of them. Using 
that number one can estimate the number of an- 
nihilations that can provide the measured UHECR 
flux, which turns out to be less than 1 per year per 
volume such as that of the Galaxy. Another pos- 
sibility is the emission of X particles from cosmic 
necklaces - a closed loop of cosmic string includ- 
ing monopoles. This particular type of topological 
defect has been extensively studied [21]. 

The other option is that the X particles themselves 
are remnants of the early Universe. Their lifetime 
should be very long, maybe longer than the age 
of the Universe |22|. They could also be a sig- 
nificant part of the cold dark matter. Being super- 
heavy, these particle would be gravitationally at- 
tracted to the Galaxy and to the local supercluster, 
where their density could well exceed the average 
density in the Universe. 

There are two main differences between bottom-up 
and top-down models of UHECR origin. The as- 
trophysical acceleration generates charged nuclei, 
while the top-down models generate mostly neu- 
trinos and 7-rays plus a relatively small number of 
protons. The energy spectrum of the cosmic rays 
that are generated in the decay of X particles is rel- 
atively flat, close to a power law spectrum of index 
a=1.5. The standard acceleration energy spectrum 
has index equal to or exceeding 2. 

Hybrid JVIodels 

There also models that are hybrid, they include el- 
ements of both groups. The most successful of 
those is the Z-burst model [23, 24] . The idea is tiiat 
somewhere in the Universe neutrinos of ultrahigh 
energy are generated. These neutrinos annihilate 
with cosmological neutrinos in our neighborhood 
and generate Zq bosons which decay and generate 
a local flux of nucleons, pions, photons and neu- 
trinos. The resonant energy for Zq production is 
4x 10^^ eV/m,XeV), where m,^ is the mass of the 
cosmological neutrinos. The higher the mass of 
the cosmological neutrinos is, the lower the res- 
onance energy requirement. In addition the cos- 
mological neutrinos are gravitationally attracted to 
concentrations of matter and their density increases 
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in our cosmological neighborhood. If the neutrino 
masses are low, of order of the mass differences de- 
rived from neutrino oscillations, the energy of the 
high energy neutrinos should increase. 

PROPAGATION OF UHECR 

Particles of energy 10^" eV can interact on almost 
any target. The most common, and better known, 
target is MBR. It fills the whole Universe and its 
number density of 400 cm^"^ is large. The in- 
teractions on the radio and infrared backgrounds 
(IRB) are also important. Let us have a look at the 
main processes that cause energy loss of nuclei and 
gamma rays. 
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Figure 3: Photoproduction cross section as a func- 
tion of the photon energy for stationary proton tar- 
gets. 



Energy Loss Processes 

The main energy loss process for protons is the 
photoproduction on astrophysical photon fields 
P7 p+riTT. The minimum center of mass energy 
for photoproduction is ^sthr — rrip+m^o ^ 1.08 
GeV. Since s = + 2(1 — cos6)Epe (where 9 is 
the angle between the two particles) one can esti- 
mate the proton threshold energy for photoproduc- 
tion on the MBR (average energy e = 6.3 x 10^^ 
eV). For cos6' = the proton threshold energy is 
Ethr = 2.3x10^" eV. Because there are head to 
head collisions and because the tail of the MBR 
energy spectrum continues to higher energy, the in- 
tersection cross section is non zero above proton 
energy of 3 x 10^^ eV. 

The photoproduction cross section is very well 
studied in accelerator experiments and is known in 
detail. Figure 3 shows the photoproduction cross 
section in the mirror system [25], as a function of 
the photon energy for stationary protons, i.e. as it 
is measured in accelerators. At threshold the most 
important process is the A+ production where the 
cross section reaches a peak exceeding 500 /ib. It is 
followed by a complicated range that includes the 
higher mass resonances and comes down to about 
100 /zb. After that one observes an increase that 
makes the photoproduction cross section parallel 
to the pp inelastic cross section. The neutron pho- 
toproduction cross section is nearly identical. 

Another important parameter is the proton inelas- 
ticity kinei, the fraction of its energy that a proton 



loses in one interaction. This quantity is energy de- 
pendent. At threshold protons lose about 18% on 
their energy. With increase of the CM energy this 
fractional energy loss increases to reach asymptot- 
ically 50%. 

The proton pair production p7 — > e+e^ is the 
same process that all charged particles suffer in nu- 
clear fields. The cross section is high, but the pro- 
ton energy loss is of order m^/nip ~ 4xlO^''E. 
Figure 4 shows the energy loss length i;oss — 
X/kinei (the ratio of the interaction length to the 
inelasticity coefficient) of protons in interactions in 
the microwave and infrared backgrounds. 

The dashed line shows the proton interaction 
length and one can see the increase of A;j„ei in the 
ratio of the interaction length to energy loss length. 
The contribution of the pair production is shown 
with a thin line. The energy loss length never ex- 
ceeds 4,000 Mpc, which is the adiabatic energy 
loss due to the expansion of the Universe for Hq 
= 75 km/s/Mpc. The dotted line shows the neutron 
decay length. Neutrons of energy less than about 
3x lO'^'^ eV always decay and only higher energy 
neutrons interact. 

The pair production process deserves more atten- 
tion since it will become important soon. Figure 5 
shows the positron spectra produced in pair pro- 
duction interactions of protons with fixed energy. 
Next to the proton energy the figure indicates the 
inelasticity coefficients in these interactions. The 
pair production cross section grows with the proton 
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Figure 4: Energy loss length of protons in interac- 
tions in MBR and IRB. The shaded area shows the 
77 interaction length. 

energy while the inelasticity coefficient decreases. 
The combination of both these parameters gener- 
ates the energy loss length that has a minimum at 
about2xl0i''eV. 

Heavier nuclei lose energy to a different process - 
photodisintegration, loss of nucleons mostly at the 
giant dipole resonance [26]. Since the relevant en- 
ergy in the nuclear frame is of order 20 MeV, the 
process starts at lower energy. The resulting nu- 
clear fragment may not be stable. It then decays 
and speeds up the energy loss of the whole nu- 
cleus. Ultra high energy heavy nuclei, where the 
energy per nucleon is higher than photoproduction, 
have also loss on photoproduction. The energy loss 
length for He nuclei in photodisintegration is as 
low as 10 Mpc at energy of 10^° eV. Heavier nuclei 
reach that distance at higher total energy. Figure 6 
shows the energy loss time of heavy nuclei - 10^^ 
s equals approximately 1 Mpc. 
UHE gamma rays also interact on the microwave 
background. The main process is 77 e+e". 
This is a resonant process and for interactions 
in the MBR the minimum interaction length is 
achieved at 10^^ eV. The interaction length in 
MBR decreases at higher 7-ray energy and would 
be about a 50 Mpc at 10^" eV if not for the ra- 
dio background. The radio background does exist 
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Figure 5: Energy distribution of positrons gener- 
ated in pair production interactions of protons of 
fixed energy in the MBR. 




Figure 6: Energy loss length of heavy nuclei (He, 
O, Si, Fe from left to right in MBR [27]. 

but its number density is not well known. Figure 4 
shows the interaction length for this process in the 
MBR and the radio background as a shaded area. 

The fate of the electrons produced in a 77 colli- 
sion depends on the strength of the magnetic fields 
in which UHE electrons lose energy very fast. The 
photon energy is than quickly downgraded and the 
77 interaction length becomes very close to the 
gamma ray energy loss length. In the case of very 
low magnetic fields (0.01 nG) the synchrotron en- 
ergy loss is low (it is proportional to ElB"^) and 
then inverse Compton scattering (with a cross sec- 
tion very similar to this of 77) and cascading is 
possible. The energy loss length of the gamma rays 
would be higher in such a case. 
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The general conclusion from tiiis analysis of the 
energy loss of protons and gamma rays in their 
propagation through the Universe is these UHE 
particles can not survive at distances of more than 
few tens of Mpc and sources of the detected cos- 
mic rays have to be located in our cosmological 
neighborhood. Every small increase of the distance 
between the source and the observer would require 
increase of the maximum energy at acceleration (or 
other production mechanism) and will affect sig- 
nificantly the energy requirement for the UHECR 
sources. 

Modification of the Proton Spectrum in 
Propagation. Numerical Derivation 
of the GZK Effect. 

Figure 7 shows the evolution of the spectrum of 
protons because of energy loss during propagation 
at different redshifts. The thick solid lines shows 
the spectrum injected in intergalactic space by the 
source, which in this exercise is 

^ = Ax E-yeME/W^-^'eV) . 

dE 

After propagation on 10 Mpc (z = 0.0025) some 
of the highest energy protons are missing. This 
trend continues with distance and at about 40 Mpc 
another trend appears - the flux of protons of en- 
ergy just below 10^'^ eV is above the injected one. 
This is the beginning of the formation of a pile-up 
in the range where the photoproduction cross sec- 
tion starts decreasing. Higher energy particles that 
are downgraded in this region lose energy less fre- 
quently and a pile-up is developed. 

The pile-up is better visible in the spectra of pro- 
tons propagated at larger distances. One should 
remark that the size of the pile-up depends very 
strongly on the shape of the injected spectrum. If it 
had a spectral index of 3 instead of 2 the size of the 
pile-up would have been barely visible as the num- 
ber of high energy particles decreased by a factor 
of 10. 

When the propagation distance exceeds redshift of 
0.4 there are no more particles of energy above 
10^^ eV independently of the maximum acceler- 
ation energy. All these particles have lost energy 
in photoproduction, pair production and adiabatic 
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Figure 7: Evolution of the cosmic ray spectrum in 
propagation through different redshifts. 

losses. From there on most of the losses are adia- 
batic. 

In order to obtain the proton spectrum created by 
homogeneously and isotropically distributed cos- 
mic ray sources filling the whole Universe one has 
to integrate a set of such (propagated) spectra in 
redshift using the cosmological evolution of the 
cosmic ray sources, which is usually assumed to be 
the same as that of the star forming regions (SFR) 
r]{z) = J?(0)(1 + z)" with n = 3, or 4 up to 
the epoch of maximum activity Zmax and then ei- 
ther constant or declining at higher redshift. High 
redshifts do not contribute anything to UHECR (z 
= 0.4 corresponds to a propagation distance of 1.6 
Gpc for Hq = 75 km/s/Mpc). After accounting for 
the increased source activity the size of the pile-ups 
has a slight increase. 

Apart from the pile-up, there is also a dip at about 
10^^ eV which is due to the energy loss on pair 
production. It is also preceded by a small pile-up 
at the transition from adiabatic to pair production 
loss. This feature was first pointed at by Berezin- 
sky&Grigorieva [28 J. 

The GZK cutoff, the pile-up and the pair produc- 
tion dip characterize the energy spectrum of extra- 
galactic protons under the assumptions of injection 
spectrum shape, cosmic ray luminosity, cosmolog- 
ical evolution and isotropic distribution of the cos- 
mic ray sources in the Universe. 



30th International Cosmic Ray Conference 



Modification of the Gamma Ray Spectra 
from top-down models. 

Because of the strong influence of the radio back- 
ground and of the cosmic magnetic fields the mod- 
ification of the spectrum of gamma rays in a top- 
down scenario is much more difficult to calculate 
exactly. There are, however, many general features 
that are common in any of the calculations. Fig- 
ure 8 shows the gamma ray spectrum emitted in a 
top-down scenario with mx = 10^^ GeV [29]. The 
spectra of 7-rays and electrons from the X decay 
chain are indicated with different line types. 
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Figure 8: Evolution of the energy spectra of elec- 
trons and 7-rays injected in a top-down scenario 
with X particle mass of 10^^ GeV. 

The typical gamma ray energy spectrum in top- 
down models is E~'^/^. We shall start the discus- 
sion from the highest energy and follow the energy 
dissipation in propagation. The highest energy 
gamma rays have not suffered significant losses. 
At sUghtly lower energy, though, the 77 cross sec- 
tion grows and the energy loss increases. One can 
see the dip at about 1Q^'^~^^ GeV which is caused 
by the radio background. The magnetic field is 
assumed sufficiently high that all electrons above 
10*^ GeV immediately lose energy in synchrotron 
radiation. The minimum ratio of the gamma-ray 
to cosmic ray flux is reached at about 10^^ eV, 
where the minimum of the interaction length in 
the MBR is, after which there is some recovery. 
There is another absorption feature from interac- 
tions on the infrared background. The gamma ray 



peak in the GeV region consists mostly of syn- 
chrotron photons. Isotropic GeV gamma rays, that 
have been measured, can be used to restrict top- 
down models in some assumptions for the mag- 
netic field strength. 

UHECR Propagation and Magnetic Fields 

The possible existence of non negligible extra- 
galactic magnetic fields would modify the propa- 
gation of the UHE cosmic rays independently of 
their nature and origin. There is little observa- 
tional data on these fields. The estimate of the aver- 
age strength of these fields in the Universe is 10~^ 
Gauss (1 nG) [30]. On the other hand /zG fields 
have been observed in clusters of galaxies, and in 
a bridge between two parts of the Coma cluster. 

Even fields with nG strength would affect the prop- 
agation of UHE cosmic rays. If UHECR are pro- 
tons or heavier nuclei they would scatter of these 
fields. This scattering would lead to deviations 
from the source direction and to an increase of 
the pathlength from the source to the observer. It 
would make the source directions less obvious and 
would create a magnetic horizon [31] for extra- 
galactic protons of energy below 10^^ eV as their 
propagation time from the source to the observer 
would start exceeding Hubble time. When the 
horizon is achieved the cosmic ray spectrum ap- 
pears flatter than it actually is. 
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Figure 9: Distance from the source that contains 
more than 50% of the emitted protons. 
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If regular magnetic fields of strength exceeding 10 
nG were present on 10 Mpc coherence length they 
would lead to significant biases in the propagated 
spectra [32] in a function of the relative geometry 
of the field, source and observer. 10^^ eV parti- 
cles would gyrate around the magnetic field lines 
and thus appear coming from a wide range of di- 
rections. The flux of such particles would be much 
higher than at injection. Only particles of energy 
above 10^° eV would be able to propagate through 
the magnetic field lines. 

Ultrahigh energy protons are also affected in prop- 
agation by the galactic magnetic fields - they scat- 
ter and acquire an angle with their direction outside 
of the Galaxy. That angle depends on the UHECR 
rigidity and direction. The largest angle should 
be when the proton has to propagate close to the 
galactic center and galactic bulge region where the 
fields, although not exactly known, are the high- 
est. Excluding the galactic center vicinity, the av- 
erage deflection angle for 10^° eV protons is be- 
tween 3.1° and 4.5° in different galactic magnetic 
field models [33]. 

Production of Secondary Particles in Prop- 
agation 

One interesting feature that can be used for testing 
of the type and distribution of UHECR sources is 
the production of secondary particles in propaga- 
tion. The energy loss of the primary protons and 
7-rays is converted to secondary gamma rays and 
neutrinos (in the case of primary nuclei). Gamma 
rays are generated in nucleon photoproduction in- 
teractions and in BH pair production processes as 
well as in 77 collisions. In the case of isotropic 
and homogeneous source distribution the gamma 
ray energy is degraded and eventually converted to 
MeV/GeV diffuse isotropic flux. The value of this 
flux could be used to restrict the amount of energy 
in UHECR [34, 29]. Some cosmologically nearby 
sources may still create a halo in the source direc- 
tion of high energy gamma rays that could be de- 
tected by the sensitive contemporary TeV gamma 
ray telescopes. 

Most interesting are the cosmogenic neutrinos, that 
were first proposed by Berezinsky & Zatsepin [35] 
and have been since calculated many times, most 
recently in Ref. [36]. Every charged pion pro- 
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Figure 10: Production of neutrinos from cosmic 
ray proton propagation on different distances from 
10 to 200 Mpc 

duced in a photoproduction interaction generates 
three neutrinos through its decay chain. 

The spectrum of cosmogenic neutrinos depends on 
the UHECR spectrum, the UHECR source distri- 
bution and very strongly on the cosmological evo- 
lution of the UHECR sources [37]. The sensitiv- 
ity to the cosmological evolution of the sources is 
very high because of the lack of energy loss (ex- 
cept for adiabatic loss) of the generated neutrinos. 
Figure 10 shows the spectra of neutrinos gener- 
ated in proton propagation at different distances. 
In the contemporary Universe muon neutrinos and 
antineutrinos peak at about 10^® eV, while electron 
neutrinos and antineutrinos show a double peaked 
spectrum. The higher energy peak, coinciding with 
that of consists mostly of electron neutrinos, 
while the lower energy one is of from neutron 
decay. If there is a high fraction of heavy nu- 
clei in the primary UHECR the would domi- 
nate the i/f^ flux since there would be many more 
neutrons from photodisintegration than photopro- 
duction interactions, which mostly protons and He 
nuclei would suffer. 

It is of some importance to note that MBR is not 
the only target for neutrino production. The sec- 
ond most important one is the isotropic infrared 
and optical background (IRB). Its number density 
is, or course, much lower, but lower energy protons 
can interact in it and even in the case of flat accel- 
eration spectra the number of interacting protons 
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Figure 1 1 : Cosmogenic + generated by pro- 
tons with a flat acceleration spectrum (7=1.0) with 
cosmological evolution and a steep one (7=1.7) 
without evolution. 

to a large extent compensates for the lower photon 
target density. 

Fig. 1 1 shows the spectra of y^+i>^ generated by 
a flat (7=1.0) and steep (7=1.7) UHECR acceler- 
ation spectra. In the case we are lucky enough to 
detect cosmogenic neutrinos with the neutrino tele- 
scopes under construction and design they could 
help a lot in limiting the models for the origin of 
the ultrahigh energy cosmic rays. 

Experimental data 

At this meeting we saw the first set of experimen- 
tal data with very good statistics observed by the 
Auger CoUaboration. Auger reported the data from 
an exposure of more than 5,165 km^.ster.yr ob- 
tained during the construction of the array. This 
exceeds the exposure of the previous biggest array, 
Agasa [38], by about a factor of three. The ex- 
posure of the HiRes detector is energy dependent 
as the fluorescent hght of higher energy showers 
can be detected from larger distances. The Auger 
air shower array consists of 1,600 water tanks of 
area 10 m^ in which the charged shower particles 
and the converted 7-rays produce Cherenkov light. 
Tanks are viewed with 3 photomultipliers. The wa- 
ter tanks have the advantage to have significant ef- 
fective area for highly incUned showers. 



Cosmic ray spectrum 

Auger reported three energy spectra obtained in 
different manner: one from the surface array nor- 
malized to the fluorescent telescopes measure- 
ment [39], another from hybrid array, i.e. from 
showers observed both by the ground array and 
by the fluorescent detectors [40], and a third one 
coming from showers arriving at zenith angles ex- 
ceeding 60° [41]. All three spectra agree with each 
other within the statistical uncertainties. There is 
only two events of energy above 10^" eV in this 
set. The presented spectra thus support the conclu- 
sion of HiRes [42, 43] that the cosmic ray spec- 
trum does not continue above IQ™ eV with the 
same ~E~^-^ spectral index as the Agasa exper- 
iment found. There is obviously a steepening of 
the spectrum which may be consistent with a GZK 
feature. 

There are good reasons to trust the spectrum mea- 
surements of the Auger collaboration. The anal- 
ysis only includes well contained showers by the 
requirement that the highest hit detector is sur- 
rounded by six active detectors. This requirement 
also guarantees that there is enough information 
for a good shower analysis. Another requirement 
is the reconstructed shower core position is inside 
the 3,000 km^ array. For this reason the exposure 
of Auger is very well known. Only events with 
reconstructed energy above 3x 10^* eV, where the 
efficiency is 100%, are included. The uncertainty 
on the energy estimate is quoted to be 22% - most 
of it due to fluorescent efflciency. 

The overall normalization of the spectrum is some- 
what lower than that of HiRes. It also seems to 
have sUghtly different shape: the dip in the spec- 
trum above 10^^ eV (seen mostly in the hybrid data 
set appears deeper and the recovery faster as shown 
in Fig. 12. The depletion of high energy showers 
is not that steep and starts a bit earlier. Note that 
the slopes presented in Fig. 12 are not the same as 
presented by Auger at the meeting [44] - my fits 
are probably not as careful as those of the group. 
As small these differences are, they currently af- 
fect the comparison with the UHECR acceleration 
and propagation models and the derivation of the 
end of the galactic cosmic ray spectrum (see the 
talk of Venya Berezinsky in this volume). The so- 
lution of these questions will not be made before 
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Figure 12: Comparison of the detailed spectral 
shapes reported by the Auger Collaboration (sur- 
face detector and hybrid events) and by HiRes 1 
and 2 in monocular mode. 

we have a very good measurement of the UHECR 
chemical composition. 

UHECR composition 

Auger also presented data on the average depth 
of shower maximum (Xmax) as a function of the 
shower energy [45]. The average X„iax is the mea- 
sure of the cosmic ray chemical composition that 
could be made by fluorescent detectors. Hybrid 
events are used in this data set because even one 
surface detector triggered in coincidence with the 
fluorescent trigger vastly improves the shower re- 
construction. Fig. 13 compares the measurements 
of HiRes, HiRes Prototype/MIA and Auger, con- 
verted to {In A) using the Sibyll2.1 hadronic in- 
teraction model. The use of other models would 
change the derived values, i.e. move the exper- 
imental points to lower {InA) values for models 
having shallower X^ax- The two sets of squares 
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Figure 13: The data of HiRes/MIA, HiRes and 
Auger are compared to the predictions of three 
models of UHECR acceleration and propagation. 

and the shaded area plot roughly the (InA) val- 
ues that come from three types of models: the 
empty squares correspond to the model of Berezin- 
sky et al [46, 47], full ones correspond to the 
model of Bahcall & Waxman [48] and the shaded 
area corresponds to the mixed primary composi- 
tion model [49, 50]. At very high energy in Fig. 13 
I have assumed a standard H + He composition in 
a ratio of 9:1. 

The model of Berezinsky et al is that all cosmic 
rays above 10^^ eV are of extragalactic origin, and 
are mostly protons. The acceleration (injection) 
spectrum is steep with a slope of 1 .7. The dip is due 
to pair production losses. The shape of the spec- 
trum below the dip is determined by the transition 
to purely adiabatic energy loss. There is no need 
for a strong cosmological evolution of the cosmic 
ray sources. The maximum amount of He and 
heavier nuclei is about 15%. Bahcall & Waxman 
assume that the dip is where the flat extragalactic 
UHECR flux (7=1) intersects the flux of galactic 
cosmic rays. A similar model was also proposed 
by Wibig & Wolfendale [5 1 ] with slightly different 
parameters. The flat acceleration spectrum models 
require strong cosmological evolution (1 + z)^^* 
to be able to fit the experimental data. Since the 
extragalactic cosmic rays are mostly H the compo- 
sition becomes very Ught after the transition. Ob- 
viously in the Berezinsky et al model there is no 
need of galactic cosmic rays of energy above 10^* 
eV, while in the other models the Galaxy has to ac- 
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celerate particles to energy higher by an order of 
magnitude. 

Mixed UHECR composition model assumes that 
these particles leave their sources with a com- 
position similar to those of GeV galactic cosmic 
rays. The composition changes in propagation and 
above 10^^'^ eV the composition also becomes 
very light. The injection spectrum is relatively flat 
(7=1.2-1.3). Since the spectrum at Earth depends 
on the primary composition as well as on the ac- 
celeration spectrum it is possible to fit the observa- 
tions in different ways. The current experimental 
data on the cosmic ray composition do not seem to 
support any of the models. 

There is, however, a problem in our understand- 
ing of the shower development in the atmosphere. 
If Auger surface detector data were analyzed us- 
ing existing shower simulations (not normalized 
to the fluorescent detector) the energy estimate 
would on average go up by 20-25% and the com- 
position would appear heavier [52]. The only 
current hadronic interaction model that predicts 
similar X^ax and high surface (muon) density is 
EPOS [53], which is not yet well studied. 

For this reason the solution of the problems of the 
highest energy cosmic rays will have to wait. The 
differences in the Auger energy estimates by the 
surface and fluorescent detectors is very similar to 
the current differences between the AGASA and 
HiRes spectra after the energy estimate of AGASA 
has come down by 10-15% with the use of the con- 
temporary hadronic interaction model [54]. The 
hybrid detection of Auger and the Telescope Ar- 
ray, as well as the theoretical work, will lead us to 
the solution in the next several years. 

One of these problems is almost solved: many 
experiments and currently Auger [55] have set 
strict limits on the fraction of gamma-rays in the 
UHECR flux. For showers above 10^^ eV the limit 
of Auger is 2%. At higher energy the Umits are 
not that strict because of limited statistics. A gen- 
eral conclusion can be drawn from these limits that 
UHECR are not the result of top-down models and 
are due to acceleration in powerful astrophysical 
objects. 



Arrival directions of UHECR 

The question then is which these objects are. The 
AGASA experiment [56] has seen clustering of the 
highest energy events, i.e. several groups of two 
events and one of three events that come from sim- 
ilar directions smaller than the angular resolution 
of the array. HiRes did not observe that type of 
clustering. Auger [57] does not confirm the clus- 
tering either, although there is still some possibility 
that some degree of clustering (2% probability for 
isotropic distribution) exists among the highest en- 
ergy events. 

So we are back to that controversial situation 
where we know that the astrophysical sources of 
UHECR have to be close to us in cosmological 
sense but we can not see them. The Auger Collab- 
oration is in the process of an intensive search for 
anisotropics and correlations with different types 
of astrophysical objects in their data sample. Since 
we expect a very strong increase in their statistics 
even during the next year we hope future data will 
help resolving this situation. 

Summary 

The cosmic ray spectrum steepens at the approach 
of 10^° eV according to the HiRes and the new 
Auger data. Is this steepening the long expected 
GZK feature or is it of a different origin? An ex- 
ample for a different origin could be the inabiUty 
of the sources to accelerate cosmic rays of energy 
higher than 5 x 10^^ eV. Solving this would require 
higher statistics than is currently available. 

The current statistics is, however, enough to estab- 
hsh the fact that l-2x 10^^ eV cosmic rays are not 
7-rays. This fact supports the acceleration bottom- 
up scenarios for cosmic ray production. At the 
higher energies more statistics is needed to set lim- 
its of a higher quality. Thus it is still possible (al- 
though not likely) that the highest energy events 
could result from top-down models. 

It is not currently obvious what the cosmic ray 
nuclear composition is in the energy range above 
10^* eV. The Auger data points at a composition 
that is somewhat heavier than the one derived by 
the HiRes experiment. Differences are not stati- 
cally significant yet. Composition measurements 



UHECR: Origin and propagation 



are hurt by the insufficient understanding of the 
shower development. Current hadronic interaction 
models should improve after comparisons with the 
Large Hadronic Colhder (LHC) data that will be- 
come available soon. 

We still can not see the UHECR sources. Studies 
of anisotropy, in addition of revealing the UHECR 
sources, are complimentary to the composition 
studies. If these particles were heavy nuclei they 
would scatter stronger in the magnetic fields and 
show smaller correlation with their sources. Pro- 
tons would more or less point at their sources. The 
identification of the UHECR sources would also 
contribute to our knowledge of magnetic fields in 
the Galaxy and the Universe. 

We are now in a period when the UHECR statistics 
is increasing very fast and will help the solution of 
the problems listed above. In addition to Auger and 
TA the work on a satellite based UHECR Observa- 
tory continues. While the previous two projects, 
EUSO [58] and OWL [59] barely exist anymore, 
the Japanese JEM/EUSO project, that is to be in- 
stalled at the Japanese module of the International 
Space Station is doing well and may be launched 
in 2013. Such an experiment would increase the 
statistics by another factor of 10 over the surface 
arrays. 

Finally, the neutrino astronomy projects are also in 
the process of fast development. Auger is look- 
ing for horizontal air showers initiated by neutri- 
nos, the construction of IceCube [60] is going ex- 
tremely well, and the radio detection of UHE neu- 
trinos is also moving forward. The possible detec- 
tion of cosmogenic neutrinos will help the solution 
of the UHECR origin when compared to the direct 
observations. 
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